Origin of complex crystal structures of elements at pressure 
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We present a unifying theory for the observed complex structures of the sp-bonded elements 
under pressure based on nearly free electron picture (NFE). In the intermediate pressure regime 
the dominant contribution to crystal structure arises from Fermi-surface Brillouin zone (FSBZ) 
interactions - structures which allow this are favoured. This simple theory explains the observed 
crystal structures, transport properties, the evolution of internal and unit cell parameters with 
pressure. We illustrate it with experimental data for these elements and ab initio calculation for Li. 
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Recent experimental high-pressure investigations of 
metallic elements have yielded surprising and intriguing 
results, with decreasing coordination and increasing crys- 
tal complexity at intermediate pressures (figQi. This be- 
havior can be reproduced by ab initio calculations based 
on the density functional plane wave pseudopotential 
method (DFPP) which represents a sufficient theoreti- 
cal understanding of the problem: there is no need for 
physics "beyond" DFPP, as is required in the /-metal 
localisation transitions. However, DFPP calculations are 
sufficiently complex that no simpler theoretical principles 
emerge. Moreover DFPP calculations require plausible 
candidate structures - in simple cases picking the "usual 
suspects" and symmetry-breaking distortions therefrom 
has worked, but when very complex structures are con- 
tenders, a more systematic approach is needed. Here, we 
demonstrate general principles which contribute to com- 
plex crystal structures, and deduce heuristics from which 
to choose candidate structures. 

DFPP gives the enthalpy of various structures, the 
lowest enthalpy structure being stable. Under pressure 
the total enthalpy comprises a band structure term (the 
eigenvalues of non-interacting electrons moving in an ef- 
fective field), an Ewald sum, pressure times volume, ex- 
change correlation and Hartree energies. The depen- 
dence on volume (fi) is as follows: for the coulom- 
bic terms (Hartree, ion- ion, ion-electron), fi" for the ex- 
change correlation (neglecting "non-linear core correc- 
tions"), Qr"^/^ for kinetic energy. Thus under pressure 
materials become more free-electron like. We note two 
aspects of DFPP. Firstly, the fact that pseudopotentials 
work indicates that repulsion between core electrons can 
be neglected. Secondly, energy minimization in DFPP 
codes is hugely improved by preconditioning^: assum- 
ing that the contribution from a plane wave basis state 
is primarily its free-electron-like kinetic energy. 

Two further effects are not explicit here: imper- 
fect screening of ionic charge as ions approach closely 
and the FSBZ energy splitting from interaction be- 
tween plane waves and ionic potentials^: AE — 



V(q.r) depends on the scattering at q=k+k' which, 
for elements, is proportional to the X-ray diffraction. If k 
falls near the Fermi level the state with increased energy 
is unoccupied, while the other is occupied. Thus FSBZ 
interaction gives a first order change in energy with crys- 
tal structure, while at other k the energy gained and lost 
cancels. FSBZ effects scale as fi"^. 

It is important to distinguish those terms which con- 
tribute most to the total energy from those which con- 
tribute to energy differences between crystal structures. 
The latter are the screened ion-ion potential and the per- 
turbation of the free electrons from FSBZ interactions. 
Given its fl^^ scaling with volume, FSBZ interaction may 
dominate at intermediate pressures with the exp[^fl^^'^) 
dependence of imperfect screening being important at the 
highest pressures. 

The central result of this paper will be the demon- 
stration that the complex elemental crystal structures at 
intermediate pressures can be simply understood using 
NFE. This theory also gives a way of picking plausible 
trial structures for total energy calculation, and explains 
a number of observed properties. To ease comparison 
with experiment, we will refer to specific interactions be- 
tween the FS and particular points in the BZ (diffraction 
peaks) rather than the equivalent|3| description in terms 
of BZ planes containing these points. 

Under pressure, group I and II metals transform from 
simple structures, bcc and fee, to complex structures. 
The striking similarity between the Cs-oC84, Rb-oC52 
and Li-c716 structures^ Q, 0| is shown in Fig|21 and in 
the diffraction patterns (FiglSJ: all are derivatives of bcc. 

Complex phase stability in Rb and CsQ, HI has been 
attributed to s — c? transfer of electrons allowing direc- 
tional bonding and hence open structures. A strong jus- 
tification for this was the now-discredited isostructural 
phase transition in Cs0 which cannot be explained in 
a free-electron picture. In a localised basis picture s ~ d 
transfer transfer certainly exists, but offers no insight into 
the nature of the structures and cannot explain the Li- 
c/16 structure(Fig[21) as there are no d-orbitals available 
to Li. DFPP calculations for c/16 Li provide an ideal 
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FIG. 1: High-pressure structural transition sequences for the 
alkah and alkaline earth elements|l,|l|l|l|ltll[ll[ll[l3, 

IT^. H^. The general pattern is similar, with low-Z materials 
having higher transition pressures. Intermediate structures 
are labelled in Pearson notation. 

testbed for the alternate FSBZ picture (FigEl: FSBZ 
effects scale strongly f2~^ with reduced volume, causing 
increases in internal parameter x, (211) band gap open- 
ing and diffraction peak intensity. At very high pressure 
this is overcome by imperfect screening of ionic charges 
(exponential in xfl^^"^). 

We have calculated the optimal value of x in c/16 
Li using DFPP.FigureQl shows variation of total energy 
with atomic position parameter x, for a range of vol- 
umes. Ewald and band structure contributions account 
for all the variation: for larger x the band structure en- 
ergy is lowered, but the Ewald energy also increases and 
ultimately becomes dominant. Thus, the band structure 
energy is responsible for complex structure stability in Li 
and [3, lis l2il I2II ■ the signature in powder diffraction is 
that one or more diffraction peaks approach fcpls^l. 

Precisely this signature exists in Cs and Rb Figur^ 
indicating similar physics. The interstitial atoms reduce 
volume per atom at the expense of splitting the the 
diffraction peak near to the Fermi level. In a real space 
picture the sequence of first increasing then decreasing 
interstitial number seems inexplicable. In the FSBZ pic- 
ture, the Fermi surface is more easily deformed at large 
Z. For Li, the PV enthalpy gained from interstitials can- 
not compensate for loss of FSBZ interaction from peak 
splitting. Rb has FSBZ interactions with three of the 
four split peaks, while Cs deforms to touch all four. 

This gives us the desired heuristic in the search for 
candidate "complex" structures for DFPP calculation: 
multiple strong diffraction peaks at the Fermi energy. It 
also suggests that a simple model of interatomic interac- 
tions in such materials should concentrate on a particular 
wavenumber in reciprocal space: a condition somewhat 
encapsulated long ago by Friedel oscillations. 

With the FSBZ picture established for alkali metals, it 
is instructive to turn to other elements. Recently, com- 




FIG. 2; Crystal structures of Li-cI16(left) Cs-III (middle) and 
Rb-III (right) viewed along the [100] axis, of the C222i unit 
cell. Interatomic distances up to 2.5A(Li), 4.7 A (Cs) and 
4.lA(Rb) are shown as solid lines. True symmetry opera- 
tors normal to the plane of projection are shown by solid 
symbols, pseudo-symmetry operators by open symbols. The 
c/16 structure can be viewed as a simple distortion of body- 
centered cubic structure, in which each atom is shifted by x 
along the bcc [111] directions in this projection we find 8- 
atom layers. For Cs-III and Rb-IIlQ identical 8-atom layers 
existll8| interspersed with 10-atom layers (black) which can 
be considered as 8-atom layers with (001) dumbbell intersti- 
tials inserted. Thus, the Cs-oC84 and Rb-oC52 structures 
are the same as Li-c/16 with density increased by interstitial 
atoms every fifth (Cs) or third (Rb) layer. 

plex Li-ty pe d istorted bcc structures have been reported 
in gallium [l7|. This extraordinary confluence of monova- 
lent and trivalent materials would appear to contradict 
the FSBZ picture. However, figOl reveals the (310)Li- 
type diffraction peak just below the Fermi level. Such 
effects have also been observed in intermediate pressure 
phases of Group IV and III-V compounds's^]. 

In Rb, Ba and Sr a new principle emerges: In real 
space BalV-type forms self-hosting "hotel" structures 
which can only be described using two interpenetrating 
crystal structures (the "host" and chains of "guests") 
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FIG. 3: Calculated diffraction patterns for bee the complex 
structures of Li, Cs, Rb and Ga. The major effect of the 
c716 distortion in Li is to throw up a 211 peak just below the 
Fermi Vector |3^ . calculated assuming a free electron sphere 
to equate to 261 = 2 sin"^ (A -^/SiV/ei^) 



[Sl l2d UM- The structures can be incommensu- 
rate, giving rise to two distinct Brillouin zones and the 
consequent additional possibihties for FSBZ interactions. 
Again their stabiUties are well reproduced by DFPP[2^. 

Figure \S[ shows their diffraction patterns - once again 
strong diffraction peaks lie near the Fermi level, and the 
FSBZ picture is dominant. Significantly, while in Ba 
and Sr there is interaction with the (201) guest, there 
is no guest reflection (with non-zero kz) near the BZ in 
Rb. Thus FSBZ interaction order the positions of adja- 
cent chains in Ba and Sr, but not in Rb: Indeed, inter- 
chain order is observed in Ba and Sr, but Rb undergoes 
a "melting" transition at low pressure. Similarly, the ra- 
tio between guest and host lattice parameters is pressure 
independent in Ba and Sr, (locked by the FS-(201)g in- 
teraction), but pressure dependent in Rb, since no FSBZ 
effects fix the guest c/a ratio. 

Other properties of the complex crystal structures in 
alkali metals relate to the FSBZ interpretation. 

The resistivity in Li[2^, rises tenfold between 40-120 
GPa, correspondingto the c/16 and other complex or- 
thorhombic phases ]^ |^ . Similar behavior is observed in 
Cs (4 GPa, the Cs-III phase) and Rb (lOGPa, the Rb- 
II and Rb-III phases) FigH shows that the FSBZ 

interaction which open pseudogaps at the Fermi level in 
Li-c/16 giving much lower electron density at the Fermi 
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FIG. 4: Top: Ab initio calculation of phonon dispersion 
relation 33] showing dynamic instability for bcc-Li at 40GPa. 
Bottom: (right) Li-DFPP energy difference between bcc and 
cI16 for various x over a range of volumes (Vb =146.4^'^, the 
volume at which cI16 is first observed). The structure is sta- 
ble between 42-50 GPa V/Vo=l-0.85. Symbols are total en- 
ergy, lines are Ewald and bandstructure only, showing that 
the Hartree and exchange-correlation contributions are negli- 
gible. The thick line indicates the optimum value of x at each 
volume, (left) Band structures along the [2x,x,x] direction for 
Vb, x=0 (lower panel) and x=0.05 (upper panel) showing the 
opening of a pseudogap due to FSBZ interactions: the surface 
touches the 24 (211) BZ planes. Elsewhere in the BZ bcc and 
cI16 band structures are similar. 



level than in fee and bcc phases. 

Superconductivity arises from coupling of electrons to 
low frequency vibrational modes; we have not performed 
detailed calculations, but complex phases tend to ex- 
hibit superconductivity through the the low frequency 
phonons associated with FSBZ effects, and phasons of 
incommensurate phases. Superconductivity is observed 
in complex structures of Li c/16 phase Cs, Ba, Sr and 
Ga. In contrast superconductivity is not detected in Rb 
up to pressures of 21 GPa 29], where the weak non-FSBZ 
interchain coupling melts the phason mode. 

Projection of our wavefunctions onto atomic orbitals 
shows increased d character with pressure, however, 
we have shown that the concept of FSBZ interactions 
provides a better simple description for the "complex" 
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FIG. 5: Calculated diffraction patterns for hotel structures 
of Ba-IV, Sr-V and Rb-IV. The tick marks show whether the 
reflection is from host (upper) or guest (lower) lattice. The 
insets show how the Brilloiun zone maintains contact with the 
Fermi surface. 



phases observed under pressure than does s — d trans- 
fer. A simple heuristic for stable phases is the existence 
of strong diffraction peaks just below the free-electron 
Fermi vector [s^. Monovalent elements achieve this by 
distortion of bcc with FSBZ at (211) peaks; divalent 
materials form host-guest structures with two Brillouin 
Zones and trivalent gallium adopts distorted-bcc struc- 
tures with FSBZ at (310) peaks. Unlike s-d transfer, the 
theory also accounts qualitatively for the chain-melting 
and pressure sensitivity of c/a in rubidium, pressure inde- 
pendence of c/a ratio in Ba and Sr, the increased super- 
conductivity transition temperature in complex phases, 
and the similarity of Li with other group- 1 elements un- 
der pressure. 

We acknowledge discussions with M.I.MacMahon, 
R.J.Nelmes and V.Degtyareva. Calculations were done 
using the VASP and PWscf packages [s^ 
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